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SPECIES DIVERSITY, ABUNDANCE AND ENVIRONMENTAL 
ASSOCIATIONS OF LOACHES (NEMACHEILIDAE AND 
BALITORIDAE) IN THE CENTRAL REGION OF THAILAND 
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ABSTRACT 

Environmental associations favorable for benthic loaches within the families Nemacheilidae 
and Balitoridae were the objective of this study that was conducted at 160 sites within five river 
systems: Southeastern, Mekong, Chao Phraya, Maeklong and Peninsula, in the central region of 
Thailand. Loaches in this study represented seven genera and 19 species; one each of Tuberos- 
chistura and Homalopteroides, two Acanthocobitis , Homaloptera and Pseudohomaloptera, four 
Nemacheilus and seven Schistura. Fewest species ( n = 3) were found in the Southeastern system 
with four species from each of the Mekong, Chao Phraya and Peninsula systems and 13 species 
in the Maeklong system. Species diversity within all systems was greatest for Schistura. Species 
and relative abundance of all fishes and of loaches varied directly and inversely, respectively, with 
river order for the two rivers in which this was examined. Abundances varied widely among spe¬ 
cies and sites with the most abundant species being Schistura aurantiaca, Schistura cf. sexcauda, 
Homalopteroides smithi and Acanthocobitis zonalternans. Species’ distributions across all five 
river systems examined by canonical correlation analysis (CCA) varied significantly with six 
environmental variables. Of greatest significance were elevation and water temparature followed 
by water velocity, ambient oxygen, pH and Silica. Several adaptations for co-existence in fast 
flowing water are identified and discussed. 
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INTRODUCTION 

Landscapes provide environmental variation that, along with biotic interactions, contribute 
to the distribution of freshwater fishes found in rivers. Fishes live not in random groupings 
but in structured associations held together by favorable environmental characteristics (Al¬ 
lan, 2004; Hayes etal ., 2006). In some cases a small number of environmental factors seem 
to exercise a strong influence on association structure while in others it is related to a wider 
range of factors (Robinson & Tonn, 1989; Edds, 1993). Often, it is assumed that structure 
and pattern in species assemblages originate from niche specialization, functional differences 
in the way species partition-limiting resources and respond to gradients in environmental and 
microclimatic conditions. Recently, this view has been challenged by neutral theory which 
considers that trophically similar species are functionally equivalent and patterns in community 
or assemblage organization can be accounted for by unusual or random processes (Hubbell, 
2001; 2005; Stokes & Archer, 2010). The current debate on these contrasting views suggests 
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that neither is sufficiently broad to account for the full range of observed patterns in natural 
communities (McGill, 2003; Volkov etal., 2007). Additionally, it appears that high diversity 
of fishes in some tropical environments may be maintained not so much by extreme speciali¬ 
zations as by frequent occurrences of restricted environmental disturbances that, in concert 
with species specific vulnerabilities, prevent species from reaching equilibrium (Connell, 
1978; Moyle & Senanayake, 1984; Watson & Balon, 1984). 

The importance of species-specific environmental factors to fish distribution has been 
demonstrated in temperate regions (Jackson et al., 2001) but by comparison these relation¬ 
ships have been largely ignored in species-rich tropical regions, particularly Southeast Asia 
(Magurran et al., 2011; Khachonpisitsak,2012) . The high diversity in Southeast Asia likely 
relates in part to the time over which specific morphologies and processes have had to evolve 
(Beaufort, 1951). On a more local scale, environmental factors such as habitat size and diver¬ 
sity are among the variables thought to be important (Gorman & Karr, 1978; Angermeier & 
Schlosser, 1989; Hugueny, 1990; Beamish etal., 2005). Chemical characteristics are known 
also to have an important influence (Fry, 1971). 

Loaches in the species-rich families Nemacheilidae and Balitoridae live mostly in 
close association with the substrate of small fast flowing rivers in Asia and Europe and 
can be abundant (Kottelat, 1990, 2001,2012, 2013). The current study expands an earlier 
study (Beamish et al., 2008) in terms of site numbers, species identity and adaptations for 
co-existence and was undertaken to better understand their occurrence and co-existence in 
rivers throughout Thailand. 


STUDY AREA 

Small rivers representing first- to fourth-order tributaries (Strahler, 1952) of larger 
rivers were sampled in the central region of Thailand between latitudes 10° 30' and 15° 30' 
North and longitudes 98° and 102° 30' East (Figure 1). Sites (n = 160) were sampled from 
rivers in five river systems following Kottelat (1989) except boundaries of the Maeklong 
and Southeastern systems following Abell et al., (2008) (Figure 1), Mekong (n = 2), Chao 
Phraya (n = 50), Southeastern (« = 5), Maeklong (« = 92) and Peninsula (n = 11) between 
October 2000 and October 2010 (See Appendix 1 for details of each site). Number of sites was 
highest in the Maeklong and Chao Phraya systems and lowest in the Mekong and Southeastern 
systems. Sites were selected as representative of landscapes in the central region of Thailand 
and provided information on environmental factors and the presence and abundance of fish 
species. Features of particular interest were elevation, canopy, river width, discharge and 
velocity in addition to water chemistry, reflecting basin geology. Length of most landscape 
sites was approximately 30 times the estimated average width to increase the probability of 
capturing all loaches and other species present (Lyons, 1992). Typically, a site was mostly 
riffles but also contained at least one small pool. Within landscapes, sites were chosen from 
road accessible locations ranging from remote, heavily forested and sparsely inhabited to 
lightly settled areas in which modest agriculture occurred to more heavily farmed or urban 
areas. Site selection presented some problems as some waterways were accessible only dur¬ 
ing a small part of the dry season, especially those at high elevations. Sites were not sampled 
closer than 150 m from the nearest bridge to avoid potential structure bias. 
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Figure 1. 


Maps of Thailand and adjacent areas (small) 
and the central region of Thailand (large). 
Heavy black broken lines on the small map 
separate borders of the six river systems within 
and surrounding Thailand; ME. Mekong; CP, 
Chao Phraya; SE, Southeastern; MK, Maek- 
long; SW, Salween; and P, Peninsula (Kotte- 
lat, 1989; Abell etal., 2008). The large map 
illustrates locations of groups of study sites, 
each indicated by a capital letter (e.g. A-X, see 
Appendix 1 for details) in the central region 
of Thailand. Dotted areas indicate reservoirs. 
Thin broken lines indicate national boundaries, 
and areas of Thailand shown in gray color in 
both maps. 
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METHODS 

Site measurements included temperature (± 0.3°C), average width (± 0.1 m), depth (± 
lcm) and mid depth velocity (± 1 cm/s), the latter three being used to calculate discharge (1/s). 
Depth of each site was the mean of 3-5 measurements made at approximately equal intervals 
across the river. Velocity was measured at the surface and adjusted to represent the mean flow 
rate (Giller & Malmqvist, 1998). Turbidity (NTU), pH (± 0.1 unit) and dissolved oxygen 
(± 0.1 mg/1) were measured with regularly calibrated meters. Ammonia (mg NH N/1) was 
measured by the salicylate method, nitrate (mg N0 3 N /l), by the cadmium reduction method, 
total iron (mg Fe/1) by the Ferro Ver method, alkalinity (as CaC'O, mg/1, pH 4.5) using the 
sulfuric acid titration method, silica (mg SiOyi) using the heteropoly method and true color 
(platinum-cobalt color units) by the platinum-cobalt method (American Public Health As¬ 
sociation, 1992). Elevation was measured by GPS (± <10m). Sites were sampled throughout 
the year although season was not included as a variable. An earlier study by Beamish et al. 
(2005) indicated seasonal changes in fish abundances and assemblage similarity from several 
streams in central Thailand to vary inversely with discharge, which was included as a vari¬ 
able in this study. 

A substrate sample was collected at each site with a hand-held acrylic corer (5 cm inner 
diameter) to a depth of 10 ± 3 cm. Particles on the substrate surface larger than the diameter 
of corer were removed and included in the sample with adjustment for their size relative to 
total sample weight. Samples were air-dried and sieved to determine particle size distribution 
by weight. Six size categories were adopted from the Wentworth scale (Giller & Malmqvist, 
1998),> 150 mm (boulder to large cobble), 150-60.1 mm (large cobble to large pebble), 60-5.1 
mm (large pebble to coarse gravel), 5-3.1 mm (medium to fine gravel), <0.5 mm (medium 
sand to silt) and the mean particle size was calculated. Substrate for each site was coded into 
one of six categories based on mean particle size with 1 being the smallest and 6, the largest. 
Substrate at a few sites was solid or almost solid bedrock and these were coded as 7. 

Fishes were captured with a Smith-Root (Vancouver, WA, USA), model 15D backpack 
electrofisher with variable output voltage (100-1100 V), pulse width (1-120 Hz) and frequency 
(100 ms-8 ms). The anode was fitted with a 280 mm diameter ring. Output voltage was varied 
inversely with water conductivity and, for the sites in this study, was mostly between 200 
and 600 V, in combination with a wave width of 60 Hz and frequencies of 1-4 ms. Settings 
were made, based on experience, to reduce damage to fish. Periodically small numbers of 
captured fishes were returned to the laboratory and held in aquaria where they exhibited no 
deleterious effects from electrofishing over periods of up to two years. 

Seines with 3 mm mesh size were installed across the upper and lower limits of a site 
and their groundlines massed with rocks to reduce the probability of emigration from or im¬ 
migration into the sample area. A site was electrofished by moving in a zigzag pattern from 
one retaining net to the other, beginning downstream or upstream based on visibility, water 
depth and velocity. Relative capture efficiency between upstream and downstream direction 
of electrofishing was compared within several larger sites and not found to differ significantly 
(ANCOVA, p < 0.05). We captured few fishes of any species < 20 mm in total length which 
was assumed to be the lower limit of vulnerability to the sampling procedure. 

Fish were anaesthetized in methane tricaine sulfonate (150 mg/1) after capture. Those 
that could confidently be identified were enumerated and, after recovery, released downstream 
from the site. When fish could not be identified in the field, a few specimens were killed by 
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an overdose of anesthetic and preserved first in 10% formalin for 7 days and then in 70% 
ethanol for subsequent identification and permanent preservation. Fishes were identified from 
a number of sources including: Smith, 1945; Kottelat, 1984,1988,1989,1990,2000,2001; 
Plongsesthee etal., 2011,2013; Rainboth, 1996; Randall & Page, 2012, 2015; Roberts, 
1982; and others. One species collected in this study from the Thong Pha Phum region of 
Kanchanaburi Province in the Maeklong system was likely misidentified as Homalopteroides 
smithi rather than by its probable correct identity, Homalopteroides modestus. Confirmation 
was not possible as specimens were not preserved for validation. Thus, for this report any 
individuals of the apparently scarce H. modestus from rivers in Thong Pha Phum were grouped 
under H. smithi. The classification system of Nelson (2006) was followed along with names 
given in Catalog of Fishes (Eschmeyer, 2015). 

Total relative abundances of each and all species within a site were calculated by the 
maximum likelihood technique based on three to four passes at each site (Carle & Straub, 
1978). When species numbers were small and not amenable to this technique, a conversion 
factor consisting of the total abundance estimate divided by total number of fish caught was 
applied to adjust numbers of each species captured. 

Associations among fish species, their abundance arithmetically adjusted to a common 
river site area of 100 m 2 , and physicochemical variables were considered using canonical cor¬ 
respondence analysis (CCA; Program CANOCO Version 4.55 October 2006 - written by Cajo 
J. F. Ter Braak (C) 1988-2006. Biometric-quantitative methods in the life and earth sciences. 
Plant Research International, Wageningen University and Research Centre Box 100, 6700 
AC. Wageningen, Netherlands. CCA was employed to identify all significant environmental 
characteristics for species and environmental variables. These relationships were described in 
an ordination diagram in which species are represented by coded numbers and physicochemi¬ 
cal variables by vectors. Statistical significance of the Monte Carlo test (1,000 permutations) 
on the relationship between species, their abundance and environmental factors was accepted 
at p < 0.05. Data were not transformed for CCA. 


RESULTS 

In the Maeklong River system, 13 loach species (Nemacheilidae, 4 species; Balitoridae, 
9 species) were found along with 71 other species from 92 sites. Longitudinal distribution of 
species in each of two rivers in this system, Khayeng (35 sites) and Pracham Mai (19 sites), 
showed numbers of loach species and those of all other fishes to be highest among the first 
order dendritic headwater streams and to decline as river order increased (Table 1). In contrast, 
relative abundance of loaches and of all other species, as number of fish/100 m 2 , increased 
directly with river order (Table 1). In each river, nine species of loaches were captured, with 
seven common to both waterways, Nemacheilidae: Acanthocobitis zonalternans , Nemacheilus 
pallidus, Schistura aurantiaca, Schistura cf. mahnerti, Schistura cf. sexcauda and Tubero- 
schistura baenzigeri ; Balitoridae: Homalopteroides smithi, (Table 2). Acanthocobitis botia 
and Schistura balteata were found only in the Pracham Mai River and Homaloptera bilineata 
(Balitoridae) and Pseudohomaloptera sexmaculata (Balitoridae) only in the Khayeng River. 
Comparative relative abundances were similar for all species present in both rivers except 
S. cf. sexcauda that was considerably the more abundant in the Khayeng River. The least 
abundant species found in both rivers were S. cf. mahnerti and T. baenzigeri. There was a 
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Table 1. Mean (± standard deviation) number of species and relative abundance (number/100 m 2 ) for 
all fish species and for loaches by river order for the Khayeng and Pracham Mai Rivers in 
the Maeklong system. 



Khayeng River (35 sites) 
Abundance Mean (± SD) 

Pracham Mai River (19 sites) 
Abundance Mean (± SD) 


Species 

Individual fish 

Species 

Individual fish 

All fish spp. 

Order 1 

10.2 (±5.0) 

262 (±203) 

8.6 (±6.8) 

273 (±98) 

Order 2 

16.5 (±3.5) 

257 (±255) 

16.0 (±0) 

115 (±261) 

Order 3 

15.9 (±4.9) 

160 (±66) 

15.3 (±1.2) 

246 (±50) 

Order 4 

23.5 (±14.1) 

117 (±80) 

26.0 (±0) 

50 (±50) 

Loach spp. 

Order 1 

1.6 (±1.5) 

15 (±19) 

2.2 (± 1.5) 

35 (±23) 

Order 2 

3.0 (±1.4) 

16 (±13) 

4.5 (±0.7) 

29 (±22) 

Order 3 

3.7 (±1.9) 

44 (±36) 

5.0 (±1.0) 

30 (±3) 

Order 4 

5.3 (±2.2) 

15 (±8) 

5.0 (±0) 

6 (± 1) 


Table 2. Range and mean relative abundance of loaches by species in the Khayeng and Pracham Mai 
Rivers in the Maeklong River system in western Thailand. SD indicates standard deviation 
of mean values. 


Species 


Abundance, number/100 m 2 


Khayeng 

(35 sites) 

Pracham Mai (19 sites) 

Range 

Mean ± SD 

Range 

Mean ± SD 

Acanthocobitis botia 

0 

0 

0-10.0 

1.2±2.8 

Acanthocobitis zonalternans 

0-26.5 

5.0±7.4 

0-22.5 

4.4±7.7 

Nemacheilus pallidus 

0-6.0 

0.5±1.2 

0-23.0 

1.8±5.3 

Schistura aurantiaca 

0-96.5 

9.0±19.4 

0-46.6 

7.6±12.7 

Schistura balteata 

0 

0 

0-2.5 

0.3 ±0.8 

Schistura cf. mahnerti 

0-13.1 

0.8±2.3 

0-2.9 

0.2 ±0.7 

Schistura cf. sexcauda 

0-119.0 

14.8±23.9 

0-18.7 

4.0 ±5.6 

Tuberoschistura baenzigeri 

0-2.0 

0.1±0.4 

0-0.9 

0.1±0.2 

Homaloptera bilineata 

0-1.5 

< 0.1±0.3 

0 

0 

Homalopteroides smithi 

0-28.5 

5.2±7.1 

0-26.4 

3.3 ±6.9 

Pseudohomaloptera sexmaculata 

0-29.4 

4.9±8.0 

0 

0 
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suggestion that A. zonalternans, S. aurantiaca and T. baezigneri occurrred in larger numbers 
at the lower order sites whereas H. smithi, S. cf. sexcauda, P. sexmaculata and, perhaps N. 
pallidus were the more abundant at the higher order sites. However, abundance by species 
varied considerably among orders, blurring patterns, if any occurred. Physical and chemical 
parameters were quite variable among sites but with similar mean values except for alkalinity 
which was the higher in the Khayeng River and turbidity and color which were higher in the 
Pracham Mai River (Table 3). 

Environmental conditions among the five river systems varied considerably, with sites 
in the western region (the Maeklong system) generally a few degrees cooler and higher in pH 
and alkalinity than in other systems (Table 4). Water in the central region was, on average, 
highest in ammonia, nitrate, total iron, color and turbidity. Physical habitat characteristics 
including canopy, discharge, water velocity and, to a lesser extent, substrate were relatively 
uniform among sites. Elevation range was greatest and least in Maeklong and Peninsula sys¬ 
tems, respectively, and similar among sites in Southeastern, the Mekong and Chao Phraya 
systems. Site values for elevation, width, depth, color and turbidity were unimodal in dis¬ 
tribution, mostly favoring lower values. Temperature was also unimodal in distribution but 
favored relatively high values. Chemical variables also were unimodal among sites. Modal 
maxima for ammonia, nitrate, iron and alkalinity occurred at comparatively low values, those 
for silica and pH at intermediate values and those for dissolved oxygen at high concentrations, 
relative to air saturation. 


Table 3. Range, mean and standard deviation (SD) of physicochemical factors for 35 sites in the 
Khayeng River and 19 sites in the Pracham Mai River, both tributaries to the Maeklong River 
in western Thailand. 



Khayeng 

Pracham Mai 

Range 

Mean ± SD 

Range 

Mean ± SD 

Elevation (m) 

165-308 

222±59 

169-853 

318±212 

Width (m) 

0.7-10.3 

5.8±2.5 

1.5-5.5 

8.0±6.6 

Depth (cm) 

5-74 

27±17 

10-49 

23±13 

Velocity (cm/s) 

5-83 

32±22 

8-55 

36± 16 

Canopy (%) 

0-90 

30±27 

0-90 

41±31 

Substrate (Coded values) 

2-7 

5±1 

2-7 

5±2 

Temperature (°C) 

17.3-28.0 

24.5±2.3 

21.8-27.3 

24.4±1.6 

Turbidity (NTU) 

2-56 

9± 12 

2-800 

89±251 

Color (CU) 

0-99 

28±44 

0-550 

70± 171 

pH 

6.5-8.3 

7.6±0.5 

6.3-8.2 

7.6±0.5 

Oxygen (mg/1) 

4.5-8.7 

7.5 ±0.9 

C/1 

SO 

1 

00 

to 

7.9±0.8 

Ammonia (mg NH.N/l) 

0-0.05 

0 

<0.01-0.08 

<0.01±<0.01 

Nitrate (mg No,N/l) 

0-4.3 

1.3±1.3 

0-1.7 

0.6±0.6 

Iron (mg Fe/1) 

0.1-0.8 

0.3±0.2 

0.0-0.79 

0.2±0.2 

Silica (mg SiO,/l) 

6.5-31.8 

19.9±7.7 

6.6-41.6 

20.0±10.6 

Akalinity (mg CaC0 3 /l) 

18-381 

135±99 

10-200 

56±53 







Table 4. Range, mean and standard deviation (SD) of physicochemical factors for 160 sites in five river systems in the central region of Thailand. 
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Loaches by species and relative abundances differed among sites and rivers. Of the 19 
species (Nemacheilidae, 13 species; Balitoridae, 6 species) captured in the central region of 
Thailand, the Maeklong system provided environmental conditions suitable for the largest 
number of species (n = 13, Tables 5-6; Figure 2). Fewest species (n = 3) were found in the 
Southeastern system. Four species were collected from each of the Mekong, Chao Phraya 
and Peninsula systems. Species diversity within all systems was greatest for Schistura with 
seven species, S. aurantiaca, S. balteata, S. crocotula , S. kohchangensis , S. cf. mahnerti, S. 
cf. sexcauda and Schistura sp. (an undescribed species). Nemacheilus contributed four spe¬ 
cies, IV. hinotatus, N. masyae, N. pallidus and N. platiceps. Acanthocobitis, Homaloptera and 
Pseudohomaloptera each contributed two species: A. botia and A. zonalternans; H. bilineata 
and H. confuzina; P. leonardi and P. sexmaculata, respectively and one each from Homalop- 
teroides, H. smithi and Tuberoschistura, T. baenzigeri. Schistura diversity was highest in the 
Maeklong system with five species as well as eight other loaches: Acanthocobitis botia, A. 
zonalternans, Homaloptera bilineata, Homalopteroides smithi, Nemacheilus binotatus, N. 
pallidus, N. placticeps and Tuberoschistura baenzigeri. Abundances varied widely among 
species, sites and river systems (Table 5) with means greatest for four species, S. aurantiaca 
(4.5/100 m 2 across all 160 sites), S. cf. sexcauda, 3.9/100 m 2 , H. smithi, 3.8/100 m 2 , and A. 
zonalternans, 3.1/100 m 2 , respectively (Table 6). 



Figure 2. Photographs of selected species: A, Schistura crocotula (total length [TL] 41 mm. Khanan 
River, Prachuap Khiri Khan Province); B. 5c/? istura cf. mahnerti (TLapproximately 68 mm, 
Khanan River, Prachuap Khiri Khan Province); C, Schistura kohchangensis (TL approxi¬ 
mately 70 mm, Pong Nam Ron River, Chantaburi Province); D, Nemacheilus masyae (TL 53 
mm, Rayong River, Rayong Province); E, Pseudohomaloptera leonardi (TL 78 mm, Pong 
Nam Ron River, Chantaburi Province); F, dorsal view of Homalopteroides smithi (TL ap¬ 
proximately 45 mm, Klong Paiboon; Chantaburi River, Chantaburi Province). Photographs 
of H. smithi and P. leonardi by R. B. Beamish; others by R. Plongsesthee. 
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Table 5. Range of loach abundance/100 m 2 by loach species within each of five river systems in the 
central region of Thailand. Numbers in parentheses represent number of sites sampled within 
each river system. ID identifies species in Figure 4. 


Range of Abundance, numbers/100 m 2 


Species 

ID 

Southeastern 

(N=5) 

Mekong 

(N=2) 

Chao Phraya Peninsula Maeklong 
(N=50) (N=ll) (N=92) 

Acanthocobitis botia 

1 





0-10.0 

Acanthocobitis zonalternans 

2 




0-16.6 

0-49.3 

Nemacheilus binotatus 

3 





0-11.6 

Nemacheilus masyae 

4 



0-4.0 



Nemacheilus pallidus 

5 





0-23.0 

Nemacheilus platiceps 

6 



0-6.2 



Schistura aurantiaca 

7 





0-96.5 

Schistura balteata 

8 





0-16.3 

Schistura crocotula 

9 




0-17.0 


Schistura kohchangensis 

10 

0-46.2 

7.3-10.7 

0-66.7 



Schistura cf. mahnerti 

11 




0-21.5 

0-49.0 

Schistura cf. sexcauda 

12 





0-119.0 

Schistura sp. 

13 





0-1.5 

Tuberoschistura baenzigeri 

14 





0-6.6 

Homaloptera bilineata 

15 





0-1.5 

Homaloptera confuzona 

16 

0-36.3 

0-57.7 




Homalopteroides smithi 

17 

1.2-46.3 

51.1-80.3 

0-3.1 

0-38.7 

0-28.5 

Pseudohomaloptera leonardi 

18 


0-62.5 




Pseudohomaloptera sexmaculata 

19 





0-46.9 


Table 6. Incidence of occurrence (%) and abundance as mean numbers/100 m 2 (SD) for each species 
for all sites (160 sites) and for those sites where one or more loach species was present (101 
sites). River systems are indicated by letters; ME. Mekong; CP, Chao Phraya; SE, Southeastern; 
MK, Maeklong; and P. Peninsula. 


Incidence Abundance, Number/100 m 2 


Species 

River 

systems 

of occurrence (%) 

Mean (SD) 

(N=160) 

(N=101) 

(N=160) 

(N=101) 

Acanthocobitis botia 

MK 

6 

10 

0.2(1.1) 

0.3 (1.3) 

Acanthocobitis zonalternans 

P, MK 

29 

46 

3.1 (7.9) 

4.5 (9.2) 

Nemacheilus binotatus 

MK 

1 

2 

0.1 (0.9) 

0.1 (1.1) 

Nemacheilus masyae 

CP 

3 

5 

<0.1 (0.3) 

0.1 (0.4) 

Nemacheilus pallidus 

MK 

10 

16 

0.3 (2.1) 

0.6 (2.4) 

Nemacheilus platiceps 

CP 

1 

2 

0.1 (0.5) 

0.1 (0.6) 

Schistura aurantiaca 

MK 

26 

41 

4.5(12.8) 

6.9(14.6) 

Schistura balteata 

MK 

3 

5 

0.2 (1.4) 

0.2 (1.6) 

Schistura crocotula 

P 

2 

3 

0.2 (1.6) 

0.3 (1.8) 
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Table 6 (continued). 


Species 

River 

Incidence 
of occurrence (%) 

Abundance. Number/100 m 2 
Mean (SD) 


systems 

(N=160) 

(N=101) 

(N=160) 

(N=101) 

Schistura kohchangensis 

SE. ME, CP 

6 

10 

1.1 (6,7) 

1.7 (7.8) 

Schistura cf. rnahnerti 

P.MK 

14 

22 

1.1 (4.8) 

2.2 (5.6) 

Schistura cf. sexcauda 

MK 

21 

33 

3.9(12.9) 

5.7(15.2) 

Schistura sp. 

MK 

1 

2 

<0.1 (0.1) 

<0.1 (0.1) 

Tuberoschistura baenzigeri 

MK 

4 

6 

0.1 (0.6) 

0.1 (0.7) 

Homaloptera bilineata 

MK 

1 

2 

<0.1 (0.1) 

<0.1 (0.1) 

Homaloptera confuzona 

SE. ME 

2 

3 

0.6 (5.4) 

0.9 (6.2) 

Homalopteroides smithi 

SE, ME, CP, 
P.MK 

39 

62 

3.8(10.0) 

5.7(11.2) 

Pseudohomaloptera leonardi 

ME 

1 

2 

0.4 (4.9) 

0.6 (5.7) 

Pseudohomaloptera sexmaculata 

MK 

16 

26 

1.4 (5.6) 

3.1 (6.7) 


Incidence of Occurrence 

Species of loaches were present in 63% of the 160 sites sampled (101 sites) in the central 
region of Thailand. Schistura were present in 45% of all sites (72 of 160) and 71% of those 
sites where loaches were found (72 of 101). Incidence of occurrence for a single species was 
highest for H. smithi (39% and 62% of 160 and 101 sites, Table 6) followed by A. zonalter- 
nans, S. aurantiaca, S. cf. sexcauda, and P. sexmaculata. Occurrences were extremely low, < 
1% for N. binotatus, N. placticeps, Schistura sp., H. bilineata and P. leonardi. 

Occurrence of some species was confined to a single river system. S. aurantiaca, S. 
balteata, Schistura cf. sexcauda and Schistura sp. were present only in the Maeklong system 
and S. crocotula only in the Peninsula system. A single species was found in two or three 
river systems. S. cf. rnahnerti was present in the Peninsula and Maeklong systems and S. 
kohchangensis in three river systems. Southeastern. Mekong and Chao Phraya. 


Environmental Associations 

In all river systems, distributions and abundance of loaches related significantly with six 
environmental variables (p < 0.05 for axis 1 and 2, respectively, Monte-Carlo test with 1000 
permutations; Figure 3). Each axis explains 10% of the variation and represents 55% of the 
species-environment relationship. The first and second axes illustrate a negative gradient of 
elevation (overall average, 237 ± 124 m), temperature (24.5 ± 3.1 ”C), velocity (35 ± 16 cm/s), 
pH (7.4 ± 0.6) and silica (22.0 ±8.3 mg SiO /I), respectively. Oxygen loaded positively (7.6 
±0.8 mg/1) on the second axis. Elevation and water temperature were of greatest importance 
to loach distribution and abundance across all river systems with water velocity, dissolved 
oxygen, silica and pH significant but of lesser importance. The significance of elevation, 
water velocity, ambient oxygen and silica were shared across all river systems with those in 
the Maeklong system, although relative importance of temperature and elevation was greater 
across all systems than in the Maeklong alone. Further, the importance of alkalinity, iron and 
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Axis 2 



Figure 3. Distribution of loaches in this study in all five river systems in the central region of Thailand 
with respect to significant environmental variables, identified by canonical correspondence 
analysis (CCA). Numbers indicate loaches identified in Table 5. Species within the smaller 
and larger circles are those that associate with average values for each of the significant 
variables for the 160 sites. Length of each environmental vector indicates the strength of its 
statistical significance relative to species association. 


river width to species distribution and abundance lessened across all river systems while that 
for pH and temperature increased relative to that in the Maeklong system. 

Over all river systems, the majority of loach species shared similar environmental con¬ 
ditions. This included A. botia, A. zonalternans, N. rnasyae, N. pallidus, S. aurantiaca, S. 
cf. mahnerti, S. cf. sexcauda, S. sp., 7. baenzigeri, Homaloptera bilineata, Homalopteroides 
smithi and P. sexmaculata. In contrast, favorable environmental conditions varied for some 
species. Schistura balteata occurred most commonly where elevations were well above the 
overall mean and ambient oxygen concentrations were well below. S. balteata was abundant 
and the only species found at the highest elevation site > 800 m. This was also the only site 
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in which aquatic vegetation was abundant. N. platiceps occurred most frequently where am¬ 
bient oxygen was high and elevation low. Schistura crocotula and Schistura kohchangensis 
tended to occupy habitats of slightly lower than average elevation and elevated oxygen but 
otherwise average conditions. Several species, N. binotatus, P. leonardi and, to a lesser extent, 
Homaloptera confuzona were associated with sites having high temperature, pH, silica and 
swift water current, all above the respective overall site means of 25° C, 7.4,22.0 mg SiO,/l 
and 35 cm/s. However, species’ occurrences and abundances for these species were low to 
suggest caution in interpretation. 


DISCUSSION 

Geographical location and morphological features of rivers contribute directly to their 
physical and chemical characteristics and these largely dictate the presence and distribution of 
available plants and animals (Beaufort, 1951). Climate and season interact to impose temporal 
abiotic and biotic changes (Gobas & McCorquodale, 1992). Current distributions of fish 
within biogeological confines are thus largely related to evolved morphological, physiologi¬ 
cal and behavioral features, processes and tolerances acting in concert with the environment. 

Information on loach habitat abounds in the literature, commonly as short accounts accom¬ 
panying species descriptions (Rainboth, 1996; Ou et al., 2011, others). Often few descriptors 
are provided and impart relatively little quantitative information (e.g. ‘shallow depth’, ‘muddy 
bottoms’, ‘fast bedrock’ and ‘torrential streams’). Dissolved oxygen, pH and temperature are 
reported sometimes. A more comprehensive approach was undertaken on Malayan balitorids 
by Alfred (1969) that improved understanding of habitat associations. He assigned associa¬ 
tions on the basis of species presence or absence into either of two categories, one based on 
substrate/vegetation and the other on water velocity with subcategories within each. This was 
supported with underwater observations. Species that occurred in several subcategories he 
called ubiquitous, those in or nearly in only one subcategory he called non-ubiquitous. With 
this Alfred assigned habitat. The study by Alfred (1969) included one species also found in 
the present study, H. smithi, to which he assigned its habitat as ubiquitous. In the present and 
in an earlier study (Beamish etal ., 2008),//. smithi was closely associated with average values 
for each of the significant factors and thus ubiquitous in Alfred’s terminology. 

A recent important and comprehensive study on habitat associations for a large number 
of fishes including loaches, was undertaken by Suvamarakshaetal. (2012). The study related 
species richness and diversity with a large number of physicochemical and geomorphological 
factors for fishes in the Ping River basin in Northern Thailand. It found six physicochemical 
factors, dissolved oxygen, temperature, pH, conductivity, phosphorus and alkalinity and six 
geomorphological factors, altitude (elevation), distance from the sea, river discharge, and 
site depth and width to be significantly related with species richness and diversity. Although 
Suvarnaraksha’s study differed in many respects from the present study both recognized the 
importance of elevation, dissolved oxygen, temperature, pH, alkalinity and site depth and 
width. It is equally clear that while these studies have contributed importantly to fish habitat 
knowledge much remains to be uncovered. 

The present study was in agreement with the earlier study by Beamish et al. (2008) in 
finding fishes in all five river systems associated with site elevation, dissolved oxygen, silica 
and ambient temperature but differing in identifying water velocity and pH as significant. 
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Most loaches occurred within a limited but similar range of environmental conditions, not far 
from their respective overall mean supporting their wide distribution among the sites sampled. 

Single species occurrences of loaches were found for several species, most often H. smithi, 
but multiple species occurrences were also found, contrary to predictions by the competitive 
exclusion principle (Gilpin & Justice, 1972; McGehee & Armstrong, 1977). This assumes 
communities exist at competitive equilibrium and requires conditions of habitat stability and 
for species to be dependent on exactly the same resources. The central region of Thailand has 
distinct rainy and dry seasons annually with rivers that undergo annual floods and droughts 
typical of pulsed ecosystems, creating unstable habitats (Junk etal., 1989; Odum etal., 1995). 
Additionally, physicochemical conditions, predation and other factors can be expected to con¬ 
tinually impose their influence on the nature of competitive interactions making it likely that 
competitive equilibrium rarely occurs, if at all (Wiens, 1977; Huston, 1979). Further, Mccann 
et al. (1998) suggest competitive equilibrium need not occur where life history variations 
differ such that an advantage at one stage in a life cycle implies a disadvantage at another. 
The high diversity of freshwater fishes in general and loaches in particular in Southeast Asian 
rivers suggests that competitive exclusion may have been partially circumvented through vari¬ 
ations in habitat stability and life histories as well as adaptations for ecological sharing, thus 
minimizing the potential negative effects of direct competition and less suitable habitats. In 
their study on resource partitioning among Sri Lankan fishes, Moyle & Senanayake (1984) 
found that microhabitat overlap was low among co-occurring species and that diets of fishes 
tended to vary when not segregated by habitat. 

The importance of biotic factors in accommodating species coexistence in riffles was 
recently demonstrated for fish by Nithirojpakdee et al. (2012,2014). Regions of high water 
velocity typically accommodate relatively large benthic populations of immature insects as 
well as other fauna that constitute a food source for loaches. Nithirojpakdee et al. (2012) 
found that while many coexisting fishes in riffles feed on the same functional feeding groups 
of macrobenthic insects (Merritt & Cummins, 1996), they display temporal and spatial differ¬ 
ences in feeding schedules that seem to foster resource sharing rather than direct competition. 
Moreover, morphologies of importance to spatial separation of habitat in flowing water can 
be inferred from shape and size of the body and fins. Loaches, particularly balitorids have 
to varying degrees compromised their shape so that it is partially dorsoventrally flattened 
lowering both drag and lift and allowing benthic refuge among rocks. Species differences in 
paired fin areas suggest a relationship with current speed and benthic station holding. Pectoral 
fin areas for P. leonardi, Homalopteroides smithi, and S. kohehangensis adjusted to a total 
length of 50 mm, were 47,43 and 38 cm 2 , respectively (F.W.H. Beamish, unpublished). In the 
riffles studied by Nithirojpakdee etal. (2012), P. leonardi occurred most frequently at micro 
riffle velocities of 175 ± 43 cm/s followed by H. smithi and S. kohehangensis at 39 ± 27 and 
34 ± 19 cm/s (F.W.H. Beamish & P. Nithirojpakdee, unpublished), in broad agreement with 
pectoral fin areas. As fish increase in size, relative surface area and hydrodynamic drag in 
flowing water can be expected to decline in contrast to a disproportionate increase in paired 
and median fin area. 

Improved understanding of lateral line systems of fish in concert with large pectoral fins 
is revealing other adaptations for specific riverine habitats. The lateral line system seems to be 
important to station-holding behaviors of fish that maintain their position in association with 
rocks (Sutterlin & Waddy, 1975). Fishes differ in the number and distribution of sensory 
neuromasts in their lateral line systems and, these, in free swimming fishes, play an important 
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role in monitoring water movement over their bodies and providing muscular feedback for 
efficient swimming in turbulent flowing water (Anderson et al., 2001; Coombs, 2001). Ar¬ 
rays of different types of sensory neuromasts interact to facilitate both orientation to water 
currents as well as to discriminate between the hydrodynamic signatures of other animals 
including potential prey (Bleckmann et al., 1991; Kanter & Coombs, 2003; Montgomery 
et al., 2003; Coombs et al., 2007). Extending large pectoral fins laterally generates negative 
lift (Coombs et al., 2007) and contributes to station holding in flowing water. It seems likely 
that lateral lines and pectoral fins of some balitorids function in a similar manner and assist 
species such as S. balteata survive at high-elevation habitats where water currents are fast 
and food resources limited. 

Loaches were associated with high dissolved oxygen despite the relatively narrow range 
found in this study from approximately 70 to 100% of air saturation. Several species were 
clearly influenced by dissolved oxygen, with N.platiceps associated with the highest oxygen 
concentrations followed by S. kohchangensis and S. crocotula, with S. balteata associated 
with lower concentrations found at high elevation sites. No information appears available on 
physiological adaptations in regards oxygen uptake within the balitorids but cutaneous respi¬ 
ration seems probable. The relationship with pH was similar to that for dissolved oxygen in 
that the range of values associated with most species was narrow and within that not usually 
considered of concern to fishes. However, blood oxygen dissociation may be impaired at low 
pH (Fromm, 1980; Alibone & Fair, 1981) which may explain a general absence of loaches 
from the lowest-pH sites. Low-pH environments also have harmful effects on ion homeo¬ 
stasis in freshwater fish (Wilson et al., 1999; Hirata et al., 2003). Silica is important to the 
distribution of grazing fishes such as balitorids that feed directly or indirectly on algae some 
of which use silica in their morphology. Most loach species were found at sites where silica 
was about 20 mg SiO.,/1. At high elevation sites silica concentrations were generally below 
10 mg SiO /I in accord with the absence of balitorids except S. balteata and N. binotatus, to 
suggest their diets are comparatively free of silicaceous algae. In contrast, H. confuzona and 
P. leonardi may be quite dependent on silicaceous algae as both were found at high silica 
concentrations, mostly above 30 mg SiCk/l. 

Argument continues on the relative importance of resource competition and sharing for 
species’ coexistence. Hubbell’s (2001,2005) neutral theory suggests that trophically similar 
species are functionally equivalent and places less emphasis on species specific habitats. Niche 
theory is more favorable towards species-specific habitats that accommodate only single spe¬ 
cies, presumably through physiological, behavioral and morphological variation in regards 
the biotic and abiotic environment. Evidence is available in support of both theories, possibly 
suggesting that neither completely accounts for all assemblages or community compositions. 


CONCLUSIONS 

Loaches are clearly masters of the substrate in fast flowing water in small rivers. This 
study examined environmental associations in sites representing diverse landscapes throughout 
five river systems in the central region of Thailand. Species were most common in Maeklong 
system in Western Thailand, their numbers and overall abundance varying directly and inversely 
with river order in accord with the general pattern for fishes other than loaches. The most 
important of the significant environmental factors across the region were elevation and ambi- 
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ent temperature. Of lesser but significant importance were dissolved oxygen, water velocity, 
ambient silica and pH. Some variation in significant environmental factors was found in sites 
in the Maeklong system. Most, but not all, loaches were associated with approximately mean 
levels of the significant factors. A few species showed strong associations with high levels 
of elevation, temperature or oxygen. The ecological importance of significant environmental 
factors to loaches relates to species-specific morphological adaptations for a benthic lifestyle 
in fast water, including body shape and fin areas and, likely, maneuverability. Potential roles 
of the lateral line system, especially the cephalic pores, for avoiding benthic obstructions and 
capturing benthic prey are considered. Although not determined directly in the current study, 
it is unlikely that the loaches studied here are tolerant of habitats in which ambient oxygen 
is much below air saturation. Rather, the charming and beautiful species within this family 
belong to small, relatively pristine and rapidly flowing streams and small rivers in which they 
are major contributors to the ecosystem. 
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Kae wanted to give so much more than time would allow. Thai fishery science has lost an 
outstanding young scientist and I, a talented young friend. 
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Appendix 1. Distribution of study sites by rivers, provinces and river systems. River names were obtained 
from maps and local residents. Study sites in close proximity to each other are grouped 
and identified by the same letter symbol. 


Site number 

Symbol 

Stream/River 

Province 

River system 

1 

A 

E-pu 

Kanchanaburi 

Maeklong 

2 

A 

Tawat 

Kanchanaburi 

Maeklong 

3 

A 

Tawat 

Kanchanaburi 

Maeklong 

4 

A 

Tawat 

Kanchanaburi 

Maeklong 

5 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

6 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

7 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

8 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

9 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

10 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

11 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

12 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

13 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

14 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

15 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

16 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

17 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

18 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

19 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

20 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

21 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

22 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

23 

A 

Pracham Mai 

Kanchanaburi 

Maeklong 

24 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

25 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

26 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

27 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

28 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

29 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

30 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

31 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

32 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

33 

A 

Ban Rai 

Kanchanaburi 

Maeklong 

34 

A 

Kratenjeng 

Kanchanaburi 

Maeklong 

35 

A 

Kratenjeng 

Kanchanaburi 

Maeklong 

36 

A 

Satamid 

Kanchanaburi 

Maeklong 

37 

A 

Pilok 

Kanchanaburi 

Maeklong 

38 

B 

Khayeng 

Kanchanaburi 

Maeklong 

39 

B 

Khayeng 

Kanchanaburi 

Maeklong 

40 

B 

Khayeng 

Kanchanaburi 

Maeklong 
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Appendix 1 (continued). 


Site number 

Symbol 

Stream/River 

Province 

River system 

41 

B 

Khayeng 

Kanchanaburi 

Maeklong 

42 

B 

Khayeng 

Kanchanaburi 

Maeklong 

43 

B 

Khayeng 

Kanchanaburi 

Maeklong 

44 

B 

Khayeng 

Kanchanaburi 

Maeklong 

45 

B 

Khayeng 

Kanchanaburi 

Maeklong 

46 

B 

Khayeng 

Kanchanaburi 

Maeklong 

47 

B 

Khayeng 

Kanchanaburi 

Maeklong 

48 

B 

Khayeng 

Kanchanaburi 

Maeklong 

49 

B 

Khayeng 

Kanchanaburi 

Maeklong 

50 

B 

Khayeng 

Kanchanaburi 

Maeklong 

51 

B 

Khayeng 

Kanchanaburi 

Maeklong 

52 

B 

Khayeng 

Kanchanaburi 

Maeklong 

53 

B 

Khayeng 

Kanchanaburi 

Maeklong 

54 

B 

Khayeng 

Kanchanaburi 

Maeklong 

55 

B 

Khayeng 

Kanchanaburi 

Maeklong 

56 

B 

Khayeng 

Kanchanaburi 

Maeklong 

57 

B 

Khayeng 

Kanchanaburi 

Maeklong 

58 

B 

Khayeng 

Kanchanaburi 

Maeklong 

59 

B 

Khayeng 

Kanchanaburi 

Maeklong 

60 

B 

Khayeng 

Kanchanaburi 

Maeklong 

61 

B 

Khayeng 

Kanchanaburi 

Maeklong 

62 

B 

Khayeng 

Kanchanaburi 

Maeklong 

63 

B 

Khayeng 

Kanchanaburi 

Maeklong 

64 

B 

Khayeng 

Kanchanaburi 

Maeklong 

65 

B 

Khayeng 

Kanchanaburi 

Maeklong 

66 

B 

Khayeng 

Kanchanaburi 

Maeklong 

67 

B 

Khayeng 

Kanchanaburi 

Maeklong 

68 

B 

Khayeng 

Kanchanaburi 

Maeklong 

69 

B 

Khayeng 

Kanchanaburi 

Maeklong 

70 

B 

Khayeng 

Kanchanaburi 

Maeklong 

71 

B 

Khayeng 

Kanchanaburi 

Maeklong 

72 

B 

Khayeng 

Kanchanaburi 

Maeklong 

73 

B 

Khayeng 

Kanchanaburi 

Maeklong 

74 

B 

Khayeng 

Kanchanaburi 

Maeklong 

75 

B 

Tuam 

Kanchanaburi 

Maeklong 

76 

B 

Bang Ka Loo 

Kanchanaburi 

Maeklong 

77 

B 

Bang Ka Loo 

Kanchanaburi 

Maeklong 

78 

B 

Bang Ka Loo 

Kanchanaburi 

Maeklong 

79 

C 

Lichia 

Kanchanaburi 

Maeklong 

80 

C 

Lichia 

Kanchanaburi 

Maeklong 

81 

C 

Lichia 

Kanchanaburi 

Maeklong 
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Appendix 1 (continued). 


Site number 

Symbol 

Stream/River 

Province 

River system 

82 

C 

Fichia 

Kanchanaburi 

Maeklong 

83 

C 

Fichia 

Kanchanaburi 

Maeklong 

84 

c 

Fichia 

Kanchanaburi 

Maeklong 

85 

c 

Thi Khrong 

Kanchanaburi 

Maeklong 

86 

D 

Kreng Kravia 

Kanchanaburi 

Maeklong 

87 

E 

Phetchaburi 

Phetchaburi 

Maeklong 

88 

E 

Phetchaburi 

Phetchaburi 

Maeklong 

89 

F 

Pranburi 

Phetchaburi 

Maeklong 

90 

F 

Pranburi 

Phetchaburi 

Maeklong 

91 

F 

Pranburi 

Phetchaburi 

Maeklong 

92 

G 

Not known 

Phetchaburi 

Maeklong 

93 

H 

Not known 

Prachuap Khiri Khan 

Peninsula 

94 

H 

Not known 

Prachuap Khiri Khan 

Peninsula 

95 

H 

Not known 

Prachuap Khiri Khan 

Peninsula 

96 

H 

Yang Khwang 

Prachuap Khiri Khan 

Peninsula 

97 

H 

Khanan 

Prachuap Khiri Khan 

Peninsula 

98 

H 

Khanan 

Prachuap Khiri Khan 

Peninsula 

99 

H 

Khanan 

Prachuap Khiri Khan 

Peninsula 

100 

I 

Ban Hin Pit 

Prachuap Khiri Khan 

Peninsula 

101 

J 

Ban Chai Thale 

Prachuap Khiri Khan 

Peninsula 

102 

J 

Not known 

Prachuap Khiri Khan 

Peninsula 

103 

K 

Kariam 

Prachuap Khiri Khan 

Peninsula 

104 

F 

Not known 

Chachoengsao 

Chao Phraya 

105 

F 

Not known 

Chachoengsao 

Chao Phraya 

106 

F 

Not known 

Chachoengsao 

Chao Phraya 

107 

M 

Not known 

Chachoengsao 

Chao Phraya 

108 

M 

Nangrong 

Nakhon Nayok 

Chao Phraya 

109 

M 

Nangrong 

Nakhon Nayok 

Chao Phraya 

110 

N 

Prachangakham 

Prachin Buri 

Chao Phraya 

111 

O 

Paknam 

Chon Buri 

Chao Phraya 

112 

O 

Tributary of Bangpakong 

Chon Buri 

Chao Phraya 

113 

P 

Kongshi 

Chon Buri 

Chao Phraya 

114 

P 

Kongshi 

Chon Buri 

Chao Phraya 

115 

P 

Kongshi 

Chon Buri 

Chao Phraya 

116 

Q 

Ban Than Trang 

Chon Buri 

Chao Phraya 

117 

Q 

Ban Than Trang 

Chon Buri 

Chao Phraya 

118 

R 

Chan Ta Than 

Chon Buri 

Chao Phraya 

119 

R 

Chan Ta Than 

Chon Buri 

Chao Phraya 

120 

R 

Chan Ta Than 

Chon Buri 

Chao Phraya 

121 

R 

Phan Sadet 

Chon Buri 

Chao Phraya 

122 

R 

Not known 

Chon Buri 

Chao Phraya 
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Appendix 1 (continued). 


Site number 

Symbol 

Stream/River 

Province 

River system 

123 

R 

Not known 

Chon Buri 

Chao Phraya 

124 

S 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

125 

S 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

126 

S 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

127 

S 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

128 

S 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

129 

S 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

130 

S 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

131 

S 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

132 

s 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

133 

s 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

134 

s 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

135 

s 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

136 

s 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

137 

s 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

138 

s 

Khao Ha Yot 

Chon Buri 

Chao Phraya 

139 

T 

Surasak 

Chon Buri 

Chao Phraya 

140 

T 

Surasak 

Chon Buri 

Chao Phraya 

141 

T 

Surasak 

Chon Buri 

Chao Phraya 

142 

T 

Surasak 

Chon Buri 

Chao Phraya 

143 

T 

Surasak 

Chon Buri 

Chao Phraya 

144 

T 

Surasak 

Chon Buri 

Chao Phraya 

145 

T 

Surasak 

Chon Buri 

Chao Phraya 

146 

T 

Surasak 

Chon Buri 

Chao Phraya 

147 

T 

Surasak 

Chon Buri 

Chao Phraya 

148 

T 

Surasak 

Chon Buri 

Chao Phraya 

149 

T 

Surasak 

Chon Buri 

Chao Phraya 

150 

T 

Surasak 

Chon Buri 

Chao Phraya 

151 

T 

Surasak 

Chon Buri 

Chao Phraya 

152 

T 

Surasak 

Chon Buri 

Chao Phraya 
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Khlong Ta Khong 
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